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1 Introduction 



It is well known that the real photon has a partonic substructure, induced by virtual fluc- 
tuations into qq pairs. These fluctuations are in part non-perturbative, and so cannot be 
calculated from flrst principles (at least not without major advances in lattice-gauge the- 
ory). Only if the parton distributions are specifled by hand at some sufficiently large input 
scale Qq can the continued evolution with be described perturbatively. Here is the 
scale of the "probing" hard process. Such a partonic substructure also exists for a space- 
like photon, = —p^ > 0. Only if is in the deeply inelastic scattering region can the 
substructure be neglected; the effects here die away like a higher-twist contribution. For 
the experimentally accessible and theoretically challenging region Aq(^^ ^ P^ ^2 GeV^ 
evolution equations for the parton distributions and their boundary conditions cannot be 
derived from perturbative QCD. In this letter we propose a theoretical ansatz to prescribe 
the modiflcation of the Q^-evolution equations of the parton distributions with changing 
P^ and the input at Qo as a function of P^. 

While the parton distribution functions (pdf's) of the real photon have been studied 
in some detail, both experimentally and theoretically (for a recent survey see e.g. [Q), 
much less is known about the virtual photon. The only published data are by the PLUTO 
Collaboration [Q. However, recently the ZEUS Collaboration presented new data from 
HERA 0. The observed distribution has been constructed for events with two jets 
above 4 GeV transverse momentum. As P^ is increased, this distribution is gradually 
suppressed at small greement with theoretical expectations. These flrst results 

will be followed by more with increasing precision. Also LEP 2 should contribute in the 
future 0. 

In a previous publication we presented several parametrizations for the parton distri- 
butions of the real photon Q , with a proposed extension to virtual photons. This is one 
of the very few studies that give explicit predictions. Drees and Godbole , following an 
analysis of Borzumati and Schuler 0, have proposed a method based on simple multi- 
plicative factors relative to the parton distributions of the real photons. A similar recipe 
has been used by Aurenche and collaborators |Q. This may be useful for estimating the 
effects of mildly virtual photons in a sample of almost real photons, but does not appear 
well suited for QCD tests of the virtual-photon distributions. A detailed study is per- 
formed by Gliick, Reya and Stratmann ^j, but a main disadvantage here is that there 
exists no explicit parametrizations of the resulting distributions. The strategy adopted is 
also only one possibility. 

In this letter we extend our previous study to a few alternative approaches for the 
virtual photon, allowing a study of the limits of our current understanding. Numerical 
approximations are introduced, which permits the parametrizations of real-photon pdf's 
to be easily extended to virtual-photon dittos. Section 2 contains a brief summary on the 
real photon, the virtual one is covered in section 3 and some comparisons are shown in 
section 4. 
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2 The Real Photon 



The parton distributions of the real photon obey a set of inhomogeneous evolution equa- 
tions: 



(91nQ2 



to leading order. The first term is the one present in standard evolution equations, e.g. 
for the proton. The second term, the so-called anomalous one, comes from branchings 
7 ~^ QQ; and is unique for the photon evolution equations. 
The solution can be written as the sum of two terms, 

Q') = /r^'l^, Q'; Ql) + Pa^'^ix, Q'; Ql) ■ (2) 

The former term is a solution to the homogeneous evolution (i.e. without the second 
term in eq. (|l|)) with a non-perturbative input at Q = Qo; and the latter is a solution 
to the full inhomogeneous equation with boundary condition f2'^'^{x,Ql;Ql) = 0. One 
possible physics interpretation is to let /J'^^ correspond to 7 ^ fluctuations, where 
V = p°, u, (p, J/ip, . . ., is a set of vector mesons, and let /J'^'^ correspond to perturbative 
("anomalous") 7 ^ qq fluctuations, q = u, d, s, c and b. The discrete spectrum of 
vector mesons can be combined with the continuous (in virtuality k"^) spectrum of qq 
fluctuations, to give 

Pai^, Q') = E ^/:'^(^, Q'l Ql) + ^E K /f ^ /r^(x, Q^; , (3) 

V JV ''Qq 

where each component /'^'^ and /'>''^^ obeys a unit momentum sum rule. Although each 
component formally depends on two scales, Ql/k"^ and Q^, it is the combination 
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that sets the length of the evolution range and thus gives the full scale dependence. 

Beyond this fairly general ansatz, a number of choices has to be made. The approach 
adopted in the SaS sets [|] is described in the following. 

What is Qo? A low scale, Qq ~ 0.6 GeV, is favoured if the V states above are to 
be associated with the lowest-lying resonances only. Then one expects Qo ~ mp/2-mp. 
Furthermore, with this Qq one obtains a reasonable description of the total 7p cross 
section, and continuity, e.g. in the primordial k± spectrum 0. Against this choice there 
are worries that perturbation theory may not be valid at such low Q, or at least that 
higher-twist terms appear in addition to the standard ones. Alternatively one could 
therefore pick a larger value, Qo ~ 2 GeV, where these worries are absent. One then 
needs to include also higher resonances in the vector-meson sector, which adds some 
arbitrariness, and one can no longer compare with low-Q data. We have chosen to prepare 
sets for both these (extreme) alternatives. 

How is the direct contribution to be handled? Unlike the p, the 7 has a direct com- 
ponent where the photon acts as an unresolved probe. In the definition of F2 this adds a 
component C^, symbolically 

F2^(x, Q') = E < [P + fl] ®C^ + P ®C^ + C^ . (5) 
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Since = in leading order, and since we stay with leading-order fits, it is permissible 
to neglect this complication. Numerically, however, it makes a non-negligible difference. 
We therefore make two kinds of fits, one DIS type with = and one MS type including 



the universal part of C^pO 



How are heavy flavours, i.e. mainly charm, to be dealt with? When jet production 
is studied for real incoming photons, the standard evolution approach is reasonable, but 
with a lower cut-off Qo ~ for 7 — > cc. Moving to deep inelastic scattering, e7 — > eX, 
there is an extra kinematical constraint: W'^ = Q'^{1 — x)/x > 4m^. It is here better to 
use the "Bethe-Heitler" cross section for 7*7* —>■ cc. Therefore two kinds of output are 
provided. The parton distributions are calculated as the sum of a vector-meson part and 
an anomalous part including all five flavours, while F2 is calculated separately from the 
sum of the same vector- meson part, an anomalous part and possibly a part now only 
covering the three light flavours, and a Bethe-Heitler part for c and b. 

Should p° and u be added coherently or incoherently? In a coherent mixture, uu : 
dd = 4 : 1 at while the incoherent mixture gives 1:1. We argue for coherence at 
the short distances probed by parton distributions. This contrasts with long-distance 
processes, such as p7 pV. 

What is Aqcd? The F2 data are not good enough to allow a precise determination. 
Therefore we use a fixed value A^^) = 200 MeV, in agreement with conventional results 
for proton distributions. 

In total, four distributions are presented 0], based on fits to available data: 

• SaS ID, with Qo = 0.6 GeV and in the DIS scheme. 

• SaS IM, with Qq = 0.6 GeV and in the MS scheme. 

• SaS 2D, with Qo = 2 GeV and in the DIS scheme. 

• SaS 2M, with Qo = 2 GeV and in the MS scheme. 

The VMD distributions and the integral of the anomalous distributions are parametrized 
separately and added to give the full result; this is of importance for the following. 



3 The Virtual Photon 

The evolution equations (in Q"^) of the pdf 's of the virtual photon (and its solutions) can 
be exactly calculated in perturbative QCD for a restricted range, namely 

Ql « « g2 . (6) 

Experimentally accessible, and theoretically challenging is, however, the low-P^ range 
Aqqj) ^ P^ ^ Ql where evolution equations cannot be derived from perturbative QCD. 
We propose pdf's that are valid for all < P^ < Q^. These are arrived at as follows. We 
start from the observation that the moments of the pdf's are analytic in the P^ plane. A 



natural way to make use of this property is |]T2[ to express them in terms of a dispersion- 
integral in the (time-like) mass square fc^ of the qq fiuctuations. This links perturbative 
and non-perturbative contributions and allows the smooth limit P^ 0. The model- 
dependence enters when specifying the necessary weight functions. We choose these in 
such a way that the resulting expressions possess the correct, known behaviours for both 
P^ and the range (^). The result is 
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A definite behaviour for lias not yet been imposed. 

Tfie non-appearance of a l/P^ contribution in (|^) can be argued in two ways. First tliis 
is wliat one expects wlien applying generalized vector-meson dominance to a continuous 
(in k"^) spectrum of qq fluctuations. Second, at large an operator-product expansion 
of the pdf's in powers of 1/P holds, but the first non- vanishing higher- twist contribution 




comes from the dimension-4 gluon condensate [T^ . 

Associating the low-A;^ part of relation (|^) with the discrete set of vector mesons gives 
a generalization of eq. (^ to 

/r (X, Q\ P') = E ^ ( zrr^ 1 ' Q'; Ql) 

2 

I ra^\x,Q';e) . (8) 

In addition to the introduction of the dipole form factors, note that the lower input scale 
for the VMD states is here shifted from Ql to some Ql > Ql- This is based on a study 
of the evolution equation J^] that shows that the evolution effectively starts "later" in 
Q'^ for a virtual photon. Qq can be associated with the Pq, Pq, Peg or Pint scales to be 
introduced below. 

Equation (|]) is one possible answer, which we will use as a reference in the following. It 
depends on both Q^ and P^ in a non-trivial way, however, so that results are only obtained 
by a time-consuming numerical integration rather than as a simple parametrization. 

In order to obtain a tractable answer, one may note that the factor (fc^/(fc^ + P^))^ 
provides an effective cut-off at A; ~ P, so one possible substitution for the anomalous 
component is 

JQI k^ 

with Pq = max(QQ, P^). All the Q^ and P^ dependence is now appearing in a combination 
like s in eq. (^, i.e. parametrizations of {x,Q'^,P'^) are readily available by simple 
modifications of those for P^ = 0. This gives the approach we adopted in ref. 0]. 

Some objections can be raised against this substitution. The choice of Pq means 
that the anomalous component is independent of P for P < Qq, and that there is a 
discontinuous change in behaviour at P = Qo- Other simple expressions could have been 
adopted to solve the problem, such a^ Pq^ = Ql + P^, but this only illustrates the 
arbitrariness of the choice. 

One guiding principle could be the preservation of the momentum sum. We recall that 











\k^ + py 









(9) 



the components f^'^ and /■T'^'i integrate to unit momentum, i.e 

y: /'dx/r(x,Q^p2)=E 



V fv Wv + Py 27r ^ Mgg F + p2 J 

(10) 

is a measure of the probability for a photon to be in a V or qq state. Momentum sum 
preservation would thus suggest the introduction of a scale P^g according to 



>l k^ \k^ + P"^ J Jp^^ k^ 
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:iii 



^ In the following the Pq prescription implies also the replacement i — > + Qo/Q^ (H) so 
as to ensure sensible behaviours for — > Ql and — > Q^. 
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which gives 



- ^ g2 + p2 I (g2 + p2)(g2 + p2) I • (12) 

For > p2 ^Yi^is simphfies to Pjj = {Ql+P'^) exp(PV(<5o+-P^))- A simple recipe is then 
to use Peff as a lower cut-off for the anomalous components (and also as the expression 
for Qo in /J'^ in (g)). 

While Pes gives the same normalization of parton distributions as does eq. (H), it does 
not give the same average evolution range and therefore not the same x shape. That is, 
eq. receives contributions from components f^'°'^{x,Q'^ , k"^) with k down to Qo, and 
thus corresponds to a larger average evolution range s than a procedure with a sharp 
cut-off k > Peff- In order to reproduce the x shape better, an intermediate Pint in the 
range Qo < Pint < PeS is to be preferred. We found no simple formula that defines an 
optimal Pint, so will use Pjnt = QoPeS as a pragmatic choice. The momentum sum can 
be preserved by a simple prefactor, i.e. in total the anomalous component is changed 
according to 



/ k' 



f, p yk^ + p2j I J ■ in(gvPi 



mt) 







lilt 





(13) 



One should note that the above approaches correspond to different evolution equa- 
tions. The pdf's of the virtual photon given by our ansatz (^ obey evolution equations 
different from those of the real photon, eq. (|l|): the homogeneous term is the same but 
the inhomogeneous one is multiplied by a factor ((^^/(Q^ + P^))^. That is, the branchings 
7* — s> qq are suppressed relative to those of a real photon, in accordance with the relation 
(0). Approximately this also holds for the Pes and Pint prescriptions. The introduction of 
Po or Pq according to eq. (||) removes the (Q^/ (Q^ -|-P^))^ factor, i.e. restores evolution to 
be fully according to eq. (|l]). Hence differences between the pdf's of the real photon and 
those of the virtual photon in the Pq and Pg schemes (and also the difference between the 
latter two) arise solely from different input distributions. Evolution equations for parton 
distributions are rigorously defined only in the range (^, where {Q'^/{Q'^ + P^))^ ~ 1. 
Hence, differences of the kind Q^/ (Q^+P^) are formally legitimate, but since the evolution 
normally is started from a Qo of the same order as P, they are non-negligible numerically. 

So far we have not imposed any constraint on the P^ — > behaviours of the pdf's 
of the virtual photon. For P^ ^ power-like terms oc (P^/Q^)^ are more important 
than the logarithmic ones (oc InQ^/P^). Then calculations based on fixed-order pertur- 
bation theory, where the full P^ dependence is kept to the order considered, are more 
appropriate than ones invoking pdf's that sum leading logarithms. Indeed, for P^ 
resolved-photon contributions originating from the quark (gluon) content of the virtual 
photon become part of the 0{as) (0(a^)) corrections^ to the leading (in ag) direct-photon 
contributions. A sensible scheme is therefore arrived at by demanding the pdf's of the 
virtual photon to approach the respective parton-model expressions, which vanish like 
InQ^/P^ for the quark distributions and faster for the gluon distribution. Such a be- 
haviour is already respected by the Pq and Pq schemes. To ensure the same limiting 
behaviour also for the other two schemes we modify the Peff and Pint schemes as follows 
(recall P^ = max(Q2, P^)): 



(l - PVQ') P^n + {PVQ') P, 



2 . 

' 



^That is, the gluon distribution vanishes faster than the quark distributions for — > 
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PL (l - PVQ') PL + (PVQ') Pl^ in In(QV^i'nt) in eq. (0) ; 

PL I — ^ (l - PVQ^) PL + (^VQ^) ^0 in the A;^ integration in eq. . (14) 

Eventually one will aime at a complete description of photon-induced reactions where 
resolved-photon contributions (i.e. those involving the pdf's of the virtual photon) are 
matched with direct-photon ones by subtracting from the latter those terms that are 
already included through the pdf's of the photon. 



4 Comparisons 

The different approaches studied above may be seen as alternatives, indicating a spread 
of uncertainty caused by our limited understanding. The PeS and Pint ones are attempts 
to obtain simple numerical approximations to eq. (|^), while the Pq and Pq ones are more 
loosely related, cf. the comment above on evolution equations. 

Further uncertainties come from the assumed pdf's of the real photon and from the 
conventional scale ambiguity problems. (Is Q set by the p_[_ of a jet or by some multiple 
thereof?) While very significant in their own right, they are not studied here. 

The above prescriptions can be applied to any set of parton distributions, provided 
that the VMD and anomalous parts have been parametrized separately. For illustrative 
purposes, we use SaS ID throughout in the following. The results are similar for set 2D 
although differences are reduced in magnitude due to the larger Qq value. 

Fig. 1 compares the u quark distribution at P^ = 0, 0.25 and 1 GeV^, for = 
10 GeV^. All approaches have in common that the small-x part of the spectrum is reduced 
more than the high-x one, proportionally speaking, reflecting the larger suppression of 
components with longer evolution range s. The Pq and Pq approaches are considerably 
above the integral of eq. (§), reflecting the difference in momentum sum. Among the 
alternatives with the same momentum sum, the shorter average evolution range for Pgfj is 
reflected in more quarks at large x (and less gluons), while Pint gives a good approximation 
to eq. (^. (In principle the integral could be performed with either scale choice for Qq in 
the VMD evolution range, but differences are minor and so we choose to use Qq = Pq as 
a simple alternative.) 

ZEUS studies the QCD jet rate for p± > A GeV, i.e. for (Q^) ~ 20 GeV^, in the 
range 0.1 ^ a; < 0.75 [^. Events with x > 0.75 are assumed to be dominated by the 
direct process, which is less sensitive to the P^ variation and is therefore used to fix the 
normalization. The ZEUS results cannot be directly compared with pdf parametrizations, 
since effects of event migration and experimental cuts are non-trivial. However, to give 
some first impression, we study the quantity 



/•0.75 f r 1 Q 1 

X(P2; = 20 GeV^) = J^^ dx lj2 + ^^(^' ^')] + l^9{x, Q'; P') , 

(15) 

which is related to the QCD jet rate. The colour factor 9/4 is the standard enhancement 
of gluon interactions relative to quark ones. The variation of X(P^)/X(0) is shown in 
Fig. 2. (It can be discussed whether one should have omitted the x factor in eq. (pTSf ); 
one would then have obtained a slightly steeper fall-off.) The integral in eq. (|^) and its 
Peff and Pint approximations give a drop by a factor of about 2 between P^ = and 
P^ = 0.5 GeV^, in rough agreement with the ZEUS data, while the Pq and Pq approaches 
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do not drop at all as much. (Note the discontinuous derivative in the Pq curve at P = Qq, 
caused by the way Pq is defined. There are also smaller kinks visible in the Pq and Pq 
curves for P ^ rric = 1.3 GeV, related to the way the charm threshold is modelled.) With 
more precise data and acceptance corrections understood, it should therefore be possible 
to discriminate among the alternatives. 

5 Summary 

In this letter we have studied the extension of the parton distributions of the real photon 
to those of the virtual one. Analyticity in P^ allows us to represent the pdf's of the 
virtual photon as a dispersion integral in the mass of the qq fiuctuations. We have 
obtained an explicit solution for the pdf's. Under the assumption of a separation of 
the qq fiuctuations in a low-mass, discrete sum of vector-meson states and a high-mass, 
continuous spectrum, various constraints on the pdf's provide us with a unique result for 
P^ ^ Q^. The final expressions for pdf's of the virtual photon Q^, P^) cannot 

be given a closed form due to their non-trivial dependence on the three variables x, Q^, 
and P^. Therefore we constructed a parametrization (the Pint prescription) which allows 
the f2 (x, Q^,P^) to be very well approximated by simple modifications to the parton 
distributions of the real photon, i.e. parametrizations are only needed for the x and 
dependence. This extension of real-photon pdf's to those of a virtual photon can be 
applied to any set of parton distributions, provided that the VMD and anomalous parts 
are available separately. It also gives F2 (x,(5^,P^). 

In order to allow for a test of the model-dependence of the pdf's of the virtual photon 
we have constructed three prescriptions alternative to Pint, all provided in closed form. 
The various prescriptions correspond to variants of order P^/ in the evolution equations 
and/or boundary conditions. The differences are readily visible in the P^ dependence of 
distributions, so HERA and LEP 2 should offer the opportunity to distinguish between 
alternatives. 

A program with the SaS parametrizations modified according to the prescriptions 
studied in this paper is available on WWW under 
http : / / thep . lu. se/ tf 2/ staff"/ tor b j or n/ Isasgain^ . 
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Figure 1: The u-quark parton distribution xu{x,Q'^ = 10 GeV^)/a;cm- Top full line for 
= 0, below comparison of five alternatives for = 0.25 and 1 GeV^. Ordered roughly 
from top to bottom, dashed is Pq, dotted is Pq, dash-dotted is Pcfr, large dots is Pint and 
full is the integral in eq. (H). 
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Figure 2: The fall-off of parton distributions with virtuality P (we have chosen P as a; 
scale rather than P^, so as to better show the small-P region), normalized to the value 
at P2 = 0, J(P2)/J(0), for g2 = 20.25 GeV^. Here J, defined by eq. (H), is the colour- 
factor- weighted sum of parton distributions in the range 0.1 < x < 0.75. Curves are 
labelled as in Fig. 1. 
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